o 



c^ 



Mon. Not. R. Astron. Soc. OOP. [TK281 (2002) Printed 2 February 2009 (MN WI^ style file v2.2) 

Dynamical Feedback of Self-generated Magnetic Fields 
in Cosmic Ray Modified Shocks 

D. Caprioli^*, P. Blasi^'^ E. Amato^ and M. Vietri^ 

^Scuola Normale Superiore, Piazza dei Cavalieri 1, 56100, Pisa, Italy 



o 
o 

(N 

O ' ^Fermilab, Center for Particle Astrophysics, Batavia, IL, USA 



^ INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi, 5, 50125, Firenze, Italy 



2 February 2009 



O 

ABSTRACT 
(N 

^ ' We present a semi- analytical kinetic calculation of the process of non-linear 

\^ ', diffusive shock acceleration (NLDSA) which includes the magnetic field am- 

(N 

^b . plification due to cosmic ray induced streaming instability, the dynamical 

t"*-- ■ reaction of the amplified magnetic field and the possible effects of turbulent 

O 

QQ ' heating. The approach is specialized to parallel shock waves and the parame- 

ters we chose are the ones appropriate to forward shocks in Supernova Rem- 
nants. Our calculation allows us to show that the net effect of the amplified 
}_j ' magnetic field is to enhance the maximum momentum of accelerated partf 



cles while reducing the concavity of the spectra, with respect to the standard 
predictions of NLDSA. This is mainly due to the dynamical reaction of the 
amplified field on the shock, which noticeably reduces the modification of the 
shock precursor. The total compression factors which are obtained for param- 
eters typical of supernova remnants are Rtot ~ 7— fO, in good agreement with 
the values inferred from observations. The strength of the magnetic field pro- 
duced through excitation of streaming instability is found in good agreement 
with the values inferred for several remnants if the thickness of the X-ray rims 
are interpreted as due to severe synchrotron losses of high energy electrons. 
We also discuss the relative role of turbulent heating and magnetic dynamical 
reaction in driving the reduction of the precursor modification. 

Key words: acceleration of particles - SNR - shock waves - magnetic field 
amplification - jump conditions 
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1 INTRODUCTION 

The supernova remnant (SNR) paradigm for the origin of galactic cosmic rays heavily relies 
on the mechanism for particle acceleration being the Diffusive Shock Acceleration (DSA) at 
the shock front generated in the supernova blast. This mechanism, also known as first order 
Fermi acceleration, has been studied in great detail and a two decades long work has led 
to the development of a kinetic non-linear theory that allows us to assess the importance 
of the dynamical reaction of the accelerated particles on the shock itself. The nonlinear 
effects described by the theory turn out to be not just corrections. They rather reveal the 
profound reasons why the mechanism works in the first place as a cosmic ray accelerator: 
the large efficiency required to explain the energetics (10 — 20% of the kinetic energy of 
the blast wave going into cosmic rays) and the large magnetic fields required to explain the 
maximum energies observed in cosmic rays are the two main motivations for developing such 
a nonlinear theory, and, not surprisingly, are among the most successful predictions of the 
theory. 



T he ini tial attempt at building a nonlinear theory led to two-fiuid models ( JDrury fc Volk 



1980 



198ll ) which provided information on the dynamics and thermodynamics of the shocked 



gas and cosmic ray gas but not on the spectrum of the accelerated particles. A satisfactory 



ki netic approach, able to predict the spectru r n of a ccele rated 



by iMalkovl fll997l ): iMalkov. Diamond fc Volkl (120001 ) and iBlasil (12002 



partic 



es was later proposed 



2004 ). More recently 



a kinetic model wh ich takes into account the possibility of arbitrary diffusion coefficients 
was put forward by lAmato &: Blasil (l2005l ). The self-excit ation of unstabl e mod es leading to 



Amato fc Blasil (1200 



magnetic field amplification was then also introduced by 

Parallel to these analytical approaches, developed primarily in the assumption of quasi- 
stationarity of the accelera tion proce s s, numerical approaches following the temporal evo- 



lution were also developed (|Bell 



1987 



Jones fc Ellison 



1991 



Ellison. Mobius fc Paschmann 



1990 



Ellison. Baring fc JonesI 



1995 



1996 



Kang fc Jones 



1997 



2005 



Kang. Jones fc Gieselei 



2OO2I ). These have been of crucial importance for the description of some aspects of the phe- 
nomenology connected with the acceleration process in SNRs, especially during the Sedov- 
Taylor phase. Numerical methods for the solution of the transport equation for cosmic rays 
and of the conservation equations for the plasma in the shock region have been exten- 
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Table 1. Parameters for 5 well studied SNRs. 



SNR 



«o(km/s) B2{[iG) 02 X 10^ 



Cas A 5200 (2500) 250-390 

Kepler 5400 (4500) 210-340 

Tycho 4600 (3100) 300-530 

SN 1006 2900 (3200) 



RCW 86 



(800) 



91-110 

75-145 



32 (36) 
23 (25) 
27 (31) 
40 (42) 
14-35 (16-42) 



The value are from IVolk. Berezhko fc Ksenofontovl 1 120051) (in parentheses) and from iParizot et alj 1120061 ). In 
ord er to estimat e the n ormalized downstream magnetic pressure 02 = -B|/(87rpoi'o)' ^°'' ^^^ SNRs discussed 
by iPar izot et alJ (l200 6p we u s ed pp = O.lmp/cm^ (SN 1006) and pg = 0.5mp/cm'^ (in the other cases). 
IViilk. Berezhko fc Ksenofontovl 1120051 ). instead, provide directly 02- 



sively used to reproduce the observed multifrequency observatio ns from single SNRs (e.g 



Volk. Berezhko fc Ksenofontov 12005 



Cassam-Chenai' et aL 



20051 ) ■ 



A typical result of all approaches to nonlinear shock acceleration is that the spectra of 
accelerated particles are far from being power laws. The concave shape leads to spectra as flat 
as p^^-"^ close to the max imum momentum, co rresponding to total compression factors that 



may exceed ~ 50 — 100 (lAmato fc Blasill2005l ). In general the total compression fact or Rtot 
is found to scale with the Mach number as Rtot oc Mq (JBerezhko fc Ellisonlll999l ). Such 



levels of shock modifi cation do not compare well with so me observations, which suggest 



Rtot ~ 7 - 10 (see e.g. 



Volk. Berezhko fc Ksenofontov 



20051). 



From the phenomenological point of view this problem has been faced by invoking some 
sort of turbulent heating: part (or all) of the energy in the form of Alfven waves which are 
responsible for the scattering of charged particles is assumed to be damped on the thermal 



gas, thereby causing its heating in the upstream precursor region (jVolk fc McKenzie 



1981 



McKenzie fc Volk 



19821 ). This process, originally investigated as a possible mechanism to 



limit the magnetic field amplification, keeping the turbulent field in the linear regime (e.g. 
6B/B ^ 1), is currently called upon also in situations in which 6B/B ^ 1. 

An important piece of information has been recently added to the debate on whether 
SNRs can be the sources of galactic cosmic rays: Chandr a X-ray obs ervations of some rem- 



Vink 



2006. and references therein 



nants showed narrow filaments of non-thermal origin (see 
for a review) . If the thickness of the rims is assumed to be due to severe synchrotron losses 
limiting the lifetime of high energy electrons, then one can infer the strength of the magnetic 
field in the downstream region, which turns out to be ~ 100 times stronger than magnetic 
fields in the ISM. 

On the other hand, it has been proposed that the narrow rims may re flect the dampiii g 
scale of the magnetic turbulence rather than the loss length of electrons (jPohl et al.ll2005l ). 



4 D. Caprioli, P. Blasi, E. Amato and M. Vietri 

This interpretation is a t odds with the morpholog y of the radio emis sion, as discussed by 



Rothenflug et aLl (J200j) (see also the discussion in lMorlino et al.ll2009l ). but at present it is 
not possible to rule it out. If it turns out to be correct, then there would be no observational 
constraint on the magnetic field in the shock region. 

If on the other hand the interpretation based on strong synchrotron losses is confirmed, 
the inferred levels of magnetization can be interpreted as a result of streami ng instabil- 



ity induced 



jy cosmic rays e fficiently accelerated at a parallel shock front (jBelll Il978 ; 



Bell fc Lucek 



2001 



Bell 



20041), alth ough alternative mechanisms have also been put for- 



ward (e.g. iGiacalone fc Jokipiil 120071 ). In Tabled] we list the SNRs where evidence has been 
collected, from X-ray observations, for a strong magnetic field: uq is the shock velocity, B2 is 
the value of the magnetic field downstream of the shock as inferred from the X-ray bright- 
ness profile and 02 is the magnetic energy density immediately downstream of the forward 
shock, in units of the total kinetic pressure at upstre am infinity (02 = B'^/i^SnpouD). The 



data are from 



Parizot et al. 



Volk. Berezhko Sz Ksenofontovl (l2005l ) (numbers in parentheses) and from 
(120061 ). 



The efficient acceleration and the magnetic field amplification are the two most impres- 
sive manifestations of nonlinear diffusive accele ration at SNR shocks. Both these a spects 
have been included in a Monte Carlo scheme by IVladimirov. Ellison fc Bykovl (l2006l ). with 



the mag netic field ampl i ficatio i i described accordingly to the pheno meno logical scenario pro- 
posed by lBell fc Lucekl (l200ll ). IVladimirov. Ellison &: Bykovl (l2006l ) find that, when amplifi- 
cation is efficient, the wave pressure makes the plasma upstream of the shock less compress- 
ible and the change in the energy density of the magnetic turbulence across the subshock 
strongly affects Rsub, which in turn affects the injection efficiency and the entire process of 
cosmic ray acceleration. 



In a previous paper (jCaprioli et al.ll2008al ) we used a three-fiuid approach (gas, cosmic 



rays and Alfven waves), to show the very general nature of the magnetic dynamical feedback: 
the shock dynamics is significantly affected by the turbulence backreaction as soon as the 
magnetic pressure becomes comparable to that of the gas upstream of the subshock. In 
particular, for the magnetization levels inferred from available observations (see Tab. [1]) and 
the saturation values that are expected from different amplification mechanisms proposed in 
the literature, we found that Rtot would naturally lead to values in the range 6-10, in good 
agreement with the values currently inferred from observations. 

In the present paper we use a kinetic model for particle acceleration at a parallel shock 
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in the non-linear regime, together with the modified conservation equations in the precursor 
and at the subshock in order to describe particle acceleration, the dynamical reaction of 
accelerated particles, the generation of magnetic field through streaming instability and the 
dynamical reaction of the magnetic field on the shock, which also results in modified cosmic 
ray spectra. 

The paper is structured as follows: in Sec. [2] we write down the correct jump condi- 



tions for a paral 
st reami ng, as in 



el, modified s 



Capriqli 



eta 



l ock wh en magnetic turbulence is excited by the cosmic ray 
(l2008al ). In Sec. [3] we summarize our kinetic model based on 
that by lAmato fc Blasil (J2006l ). including the self-consistent treatment of the magnetic field 
amplification via resonant streaming instability. The latter is described in detail in Sec. IH 
In Sec. |5]we present our main results, namely the solutions for DSA with resonant stream- 
ing instability: in particular we discuss the reduced modification of the precursor and the 
consequent steepening of the spectrum near the (increased) maximum momentum. In Sec. [6] 
we investigate the combined effect of magnetic reaction and turbulent heating and we show 
that the dominant effect on the precursor modification is likely to be that of the magnetic 
reaction. Our conclusions are in Sec. [71 



2 DYNAMICS OF A MAGNETIZED COSMIC RAY MODIFIED SHOCK 

The pressure of the accelerated particles upstream of the shock surface leads to the formation 
of a shock precursor, in which the fiuid speed gradually decreases while approaching the 
shock. One can describe this effect by introducing two compression factors Rtot = uq/u^ 
and Rsuh = Wi/m2, where u is the fluid velocity and the indexes '0', '1' and '2' refer, here 
and in the following, to quantities taken at upstream inflnity (x = — oo), and immediately 
upstream {x = 0~) and downstream {x = 0"*") of the subshock, respectively. 

We consider a non-relativistic plane shock whose normal is parallel to the background 
magnetic fleld Bq. The equations deflning the jump conditions at the shock surface in the 
presence of cosmic rays and self-generated Alfven waves can be written as 

H? = , (1) 

[pu^+p + p^]l = 0, (2) 

- 2 

(3) 



1 3 7 

2 7 — 1 



up + F^ 







where p, u, p and 7 stand for density, velocity, pressure and ratio of speciflc heats of the 
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gas, Pw and F^ are the magnetic pressure and energy flux and the brackets indicate the 
difference between quantities downstream and upstream of the subshock {[X]l = X2 — Xi). 
Some subtleties of the treatment of mass, mo mentum and energy con s ervatio n in cosmic ray 
modified shocks are extensively discussed in ICaprioli. Blasi. Amatol (J2008bl ). Here we only 
want to emphasize that the contribution of the terms related to the pressure and energy 
flux of the cosmic rays {pc and F^) disappears when considering the subshock (i.e. Pc and 
Fc are both continuous across the subshock). Also the contribution to the pressure by the 
background magnetic field vanishes, Bq being parallel to u and therefore unaffected by any 
fluid compression. The subshock Rankine-Hugoniot relations are therefore not affected by 
the contribution of cosmic rays, but they do take into account the presence of magnetic 
turbulence, leading to a magnetized gaseous shock. 
In order to infer the magnetic f i eld jum p conditions, we use the approach ofC 



Scholer fc Belcher 



(JI971I ) and IVainio fc Schlickeiserl (119991 ) to describe the transmission and reflection coeffi 



Vainio fc Schlickeiser 



dents appropriate for Alfven waves at the shock surface. Following 
(I1999I ). we account for two upstream wave trains with helicities H^ = ±1, and for their 
respective counterparts downstream. It is worth stressing that, in principle, this part of 
the problem can be a complex one, because waves can be transmitted and reflected also 
within the precursor, due to the gradient in all quantities there. In general this could lead 
to isotropization of the turbulence and to the generation of a wave train propagating in the 
direction opposite to that directly excited by the cosmic ray streaming (i.e. the one with 
He = —1). However, we do not expect this effect to be very relevant because the shocks 
we are dealing with are not very modified, as we can check a posteriori. For these reasons 
we neglect turbulence isotropization in the precursor, but we include the transmission and 
reflection of Alfven waves at the subshock surface, which enter the jump conditions. 

Let 6Bf^ be a given mode, indicated by the subscript /i, of the magnetic turbulence. We 
write the corresponding velocity perturbation simply as 

6u=-H,.^^. (4) 

Neglecting the electric field contribution, which is of order v? /c^, the magnetic pressure and 
the energy flux are respectively 



^Av') 



p^ = ^\yyBA , (5) 
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2 



E 



47r 



u 



+ Hc,^Va) + 



T.JB, 



-u . 



in 



(6) 



where the sum over fi has to be intended as the sum over all the wave modes present at a 
given position. 

In obtaining these relations we explicitly used the fact that, for Alfven waves, F^ in 
Eq. [6] has two contributions: the former is the normal component of the Poynting vector 
{B X 6u + 6B X u) X 6B/Att, while the latter represents the kinetic energy flux in transverse 
velocity, p/26u'^u (with 6u given by Eq. Hj). If the turbulence manifests itself in forms other 
than resonant Alfen waves, pw and F^ may differ significantly. We allow the generated waves 
to have both helicities, hence the sum over the modes should account, in the upstream region, 
for both forward- and backward-going waves, namely 6Bij-. Each of these modes, however, 
when advected behind the subshock has to split into two waves with opposite helicities 
in order to satisfy the Maxwe ll equations at the subshock (see e.g iMcKenzie fc Westphal 



1969 



Scholer fc Belcher 



19711 ). This fact leads to four different downstream modes, which 
can be described by the introduction of proper transmission and reflection coefficients. In 
this scenario a "reflected" ("transmitted") wave is a wave with a helicity opposite (equal) 
to the one of its upstream counterpart. It is worth stressing that, the fluid being super- 
Alfvenic either upstream or downstream, all the modes are actually advected with the fluid, 
independent of their direction of propagation. 

These coefficients, T and TZ respectively, were derived by 



McKenzie fc Westphall fll969l ). 



and for each upstream helicity H^ = ±1 read: 
Rsub + VRsub Mai + He 



T 



n 



-'^sub ' V ^sub 



R 



sub 



2 Mai + yRsubHc 

- y/Rsub Mai + He 



R 



sub 



VR. 



■sub 



(7) 



(8) 



2 Mai — VRsubHc 2 

Here Ma = u/va is the Alfvenic Mach number, namely the ratio between fluid and Alfven 
speed, which is of order 100 and more for a typical supernova shock. For each sign of H^ we 
have 6B2/SB1 = T + 71 = Rsub, and hence 



Pw2 = PwlR 



sub 



(9) 



Since the subshock can be viewed as a magnetized gas shock, as we already stressed, the 



jump conditions found by 



Vainio fc Schlickeiserl (119991 ) for the pressure and temperature 



hold also in our case and can be written respectively as: 
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P2 (7 + l)Rsuh - (7 - 1) + (7 - l){Rsub - 1)A 



Pi 



7 + 1 - (7 - l)R,ub 



P2 1 



Pi Rsub 

with A defined as: 



A 



^sub ~r i [PiiiJi ^-tisub [-'^uiJi 



(10) 

(11) 

(12) 



Rsub - 1 Pi Rsub - 1 PlUl 

Using the expressions in Eqs. [7] and [8] for the transmitted and reflected Alfven waves, we 
find: 



A 



[-tisub ~ J- J h Usub 



-Pl- 



-?^rPi 



(13) 



Following IScholer fc Belched (jl97ll ) and IVainio fc Schlickeised (jl999l ). we assume that 



the two opposite-propagating waves carry magnetic fields SBi± displaced in such a way that 
6Bi_ ■ 6Bi^ = 0. This is not the most general possible configuration, but still it is indeed 
expected to be the most common, since it describes situations in which there is only one wave 
train, or the two fields are reciprocally orthogonal, or even, on average, when the relative 
phase between the wave trains is arbitrary. 

Hereafter we use quantities normahzed to the values of ram-pressure and velocity at 
upstream infinity: 



U{x) 



u{x) 



a{x) 



(E,^4) 

SttpoMq 



P{x) 



p{x) 



^x) 



Pc{x) 

PoMq 



(14) 



If the heating of the upstream plasma is due only to adiabatic compression, using the mass 

conservation p{x)u{x) = poUo, the normalized plasma pressure can be written as 

U{x)~^ 



P(x) 



7M2 



(15) 



where as usual Mq is the sonic Mach number at upstream infinity. 

Substituting Eq. [10], Eq. [13] and the above expression for P{x) in the equation for momen- 
tum conservation, the compression factors Rsub and Rtot are related through the equation 



R, 



7+1 
tot 



^^^O^sub 



7 + 1 - Rsubil - 1) 



1+Ab 
which is the same as the standard relation (see e.g. 

where 

KB = W[l+Rsub{2h-l)] , 
and we have defined 
W = ai/Pi . 



(16) 



Blasi 



20041 ) apart from the factor (l+A^), 

(17) 

(18) 
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It is clear that the net effect of the magnetic turbulence is to make the fluid less compressible: 
if VT^l, Rtot may be considerably reduced , while the pre s sure ( and temperature) jump 



increases, according to Eq. [TO] (and [TTT) . In ICaprioli et al.l (l2008al ) we showed, by means 



of purely hydrodynamical considerations, and without any assumptions on the details of 
particle acceleration and magnetic field generation, that this is very likely the case for the 
SNRs listed in Tab. [U 

As a final remark, we notice that if one naively assumed that downstream Fyj2 = 3m2Pu;25 
instead of using the appropriate transmission and reflection coefficients (which come from 
the need of satisfying Maxwell equations at the subshock), one would obtain 

A' = [{R,^, - If - 2R,^,]W < A . (19) 

Using A' rather than A leads to an incorrect estimate of the pressure jump (Eq. [TOl) . which 
in this case may even turn out to be smaller than for an unmagnetized shock (A' < for 
Rsub < 3.73). At the same time, one would have A'^ = VT [1 + Rg^b (3/7 — 2)] in Eq. [16] and 
hence find a less marked decrease of Rfot- 

3 THE KINETIC SELF-CONSISTENT SOLUTION 

In this section we describe the calculations that lead to an exact solution for the spectrum 
and the spatial distribution of the particles accelerated at a non-linear astrophysical shock, 
including the generation of Alfven waves by the same particles and the dynamical reaction 
of both cosmic rays and magnetic turbulence on the fluid. This m ethod is based o n the 



kinet ic treatment of the problem in the stationary regime proposed by Amato fc Blasi] (J2005I . 



20061 ). which also allows for an arbitrary choice of spatial and momentum dependence of the 



diffusion coefficient D{x,p). In the following we consider Bohm diffusion in the self-generated 
magnetic field, i.e. we set 

D(.,p)^lcr,m-lc-^y (20) 

where r^ is the the Larmor radius of a particle of momentum p in the local amplified magnetic 



field B{x) = ^y8'^^a{x)pQUQ. Needless to say that this form of the diffusion coefficient is 
basically only an ansatz and whether Nature provides such a diffusion coefficient is at present 
not clear. 

The transport of accelerated particles consists of both advection and diffusion. The cor- 
rect treatment of these processes in the limit of small perturbation amplitudes {6B ^ Bq), 
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also incl uding; the presen ce of both wave trains, are well known and can be found e.g. in the 
work by ISkilling) (jl975al jbl). Unfortunately, in the case of strong magnetic field amplification 
a full theory of cosmic ray transport is still missing, and even the definition of an effective 
wave velocity, va, is troublesome. A semi-analytical treatment, which is what we are inter- 
ested in, is only possible within the framework of quasi-linear theory. We adopt an Alfven 
velocity defined as 
Bo 



VAix) 



(21) 



where p{x) is the gas density in the precursor at the position x and Bq is the strength of the 
unperturbed magnetic field. In the unlikely case that the waves could keep their Alfvenic 
nature even in the strong turbulence regime, this wave velocity would in fact be well defined. 
As to the interactions between streaming particles and Alfven waves, we do not neglect 
the velocity of the scattering centers va with respect to the fluid velocity u{x). This means 
that, in principle, the compression ratios of the background fluid are different from the ones 
experienced by the scattering centers, which turn out to be 

Ssu, = ^^i^^ S,,, = "^^^^^ ^ ^^L . (22) 

U2 + VA2 U2 + VA2 ^2 + VA2 

In the recipes above we have implicitly assumed that upstream waves are produced prefer- 
entially with H^ = —1. Behind the subshock, instead, the net velocity of the two opposite- 
propagating wave trains, when calculated in the downstre am gas ref erence frame (f ^2), turns 



Bell 



197: 



out to be in the same direction as the fluid one (see also 

In a typical SNR the condition f ^ ^ m usually holds, hence the difference between 
{Ssub, Stot) and {Rsub, Rtot) is expected not to be very relevant. But if for some reason Ma is 
small enough, the compression ratios felt by the accelerated particles may be signiflcantly 
different with respect to the fluid ones, leading to a modifled spectral slope, as already 
showed by iBelll (119781 ). This is why we retain va in the calculations and check a posteriori 
that in the cases considered this correction is not important. In Sec. 15.11 we investigate the 
consequences of adopting a different prescription for the velocity of the scattering centers. 

From the kinetic point of view, cosmic rays are described by their distribution function 
in phase space /(x, p). Keeping only the isotropic part (since f{p) = f{p) + 0{u^/(?)) 
and recalling that the shock is non-relativistic, the diffusion-advection equation for a one- 
dimensional shock reads: 



[u{x) -va{x)] 



dfix,p) 
dx 



d_ 

dx 



d 

D{x,p)—f{x,p) 



, d[u{x)-VA{x)] p df(x,p) ,^^^ 
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with Q{x,p) the injection of particles in the accelerator. We assume that injection occurs 
only at the shock location (x = 0) and at momentum Pmj, involving a fraction t] of the 
particles crossing the shock, such that 

PlUi 



Q{x,p) =7] 



Airnipp: 



2 6{p - Pinj)6{x) . 
'inj 



(24) 



For the fraction t] of injected particles we adopt the recipe of iBlasi. Gabici fc Vannoni 



fl2005h : 

4 
V 



30F 



{Ssub - l)^'e-^" 



(25) 



which assumes that only particles with momentum pinj ^ '4'Pth,2 (i-e. ip times the thermal 
particles' momentum downstream) can be accelerated. This recipe fits well within our self- 
consistent approach because it only involves pth,2 and Sgub, which are both outputs of the 
calculation, rather than free parameters. Moreover, it self-limits the acceleration process, 
suppressing the injection when particle acceleration is too efficient, in which case the large 
shock modification leads to Ssub -^ 1 and 77 -^ 0. In the following, we consider values of ip 
between 3.5 and 4, correspondin g to y betwe en ~ 10~^ an d ^ 10~^. 



Blasi. Amato fc Capriolil (120071 ). 



As shown by lAmato fc Blasil (120051 . l2006l ) and then by 
a very good approximation for the solution of Eq. [23l f{x,p), is found in the form: 



f{x,p) = /i(p)exp 



Ssub - 1 q{p)uo 



where Va 



va/uo, /i = /(0,p) and q{p) 



D{x',p) 



(26) 



d log /i (P) 
dlogp 



is the spectral slope at the shock 



location. The above expression reduces to the correct distribution function in the test particle 
limit and exactly satisfies the jump conditions at the subshock, as obtained by integrating 
Eq. [23] from 0" t o 0+. 



(27) 



As shown by iBlasil (120021 ) , /i (p) can be written as 



flip) 



3S 



tot 



VPi 



,StotUp{p) -1/ 47rmpp;„^ 

where we defined 

1 



3 exp 



P dp' 3StotUp{p') 
p,„, P' StotUpip') - 1 



Up{p) = Ui- Vai - 



^x/(x,.)^^(^)-^-(^)] 



/i(p) J-oo""" ■""'"' dx 

The normalized pressure in cosmic rays is written in terms of this solution as 

A'TT- /'Pmax 



ax) 



3poMo ^P,„, 



dpp V (p) /i (p) exp 



\^, Uix')-VAix') 
Xp{x ,p) 



(28) 



(29) 



having defined 
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Ssub - 1 Uoq{p) 
and Pmax as the maximum momentum achievable by th e accelerated particles. We determine 

the latter self-consistently, using the calculations by iBlasi. Amato fc Caprioli! (120071 ) for 

cosmic ray modified shocks, assuming that the particles' maximum energy is limited by 

the acceleration time rather than by the size of the system. The time needed to accelerate 

particles up to Pmax turns out to be 



^{Pn 



3R 



tot 



Un 



dp' 



Pinj 



p' RtotUpip') 



RtotD2{p) + 



uO 
W) 



dxf{x,p') 



(31) 



Integrating by parts Eq. [28] it is possible to express Up{p) in terms of U{x) and Xp{x,p) 
alone: 

,u{x')-VA{x'y 



Upip) 



[U{x) - Va{x)Y 

dx- -. r — — exp 

Xp{x,p) 



dx 



Xp{x',p) 
with respect to x, we obtain 

e(x)A(x) [U{x) - Va{x)] , 



Then, differentiating Eq. 
d^x) 

dx 
where 



A(x) 



j^2Tdpp''v{p)^-^^f{x,p) 
r""''"dpp^v{p)f{x,p) 

^ Ptnj 



(32) 



(33) 



(34) 



Finally, we need a relation which describes how the Alfven waves are excited by the streaming 
cosmic rays and how the wave energy is transported in the precursor. Very generally, we 
can assume a{x) to be a function of ^{x) and U{x). Using the equation of momentum 
conservation between a point x in the precursor and upstream infinity, it is therefore possible 
to write a{x) as a function of U{x) only (see e.g. Eq.|42]in the next section, where a discussion 
of magnetic field amplification due to resonant streaming instability will be presented.) 

The nonlinear system defined by Eqs. [23ll34l can be solved, for a given age of the system, 
with three nested iterations. We guess a value for Pmax-, as the starting point of the outermost 
cycle. Then, we fix a value for the ratio Rsub/Rtot, and derive the corresponding Rgub and 
Rtot from Eq. [161 The equation for conservation of momentum, 

1 



U{x) + ^{x) + a{x) + P{x) = 1 + 
once it is evaluated in 0~, 



iM^ 



2 ' 



^sub -L 



R. 



tot 



+ 6 + «1 = 1 + 



1 



(35) 



(36) 



Rtot ' iMi \R,^J ^' -fM^ ' 

only involves ^i and Rgub/Rtot, since ai can be written as a function of Rsub/Rtot, as it will 

be shown in the next section (Eq. [^2l) . Eq. [36] gives the boundary condition ^i for Eq. 
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The latter is then solved recursively: the solution at step n is calculated using U{x) and 
A(x) at step n — 1: 



e 



Wi 



XI 



6exp 



dx'X^''~^\x') [U{x') - Va{x')]'^"~'^ 



Uo 



(37) 



The iteration on n is carried out until convergence is reached. The functions U{x) (and thus 
Ssub and Stot), Up{p), A(x) and fi{p) are recalculated at every step, through Eq. [35], Eq. [32l 
Eq. [3l] and Eq. [27] respectively. In this way the value of ^(x = 0) obtained by integration 
of fi{p) through Eq. [22] in general does not match the value of ^i given by Eq. [33 hence 
we restart the procedure with a different ratio Rsub/Rtot until this necessary condition is 
satisfied. Having thus found a solution, one is able to calculate the acceleration time which 
corresponds to it according to Eq. [HU Adjusting Pmax obviously allows us to find the whole 
self-consistent solution for a fixed age of the SNR. 



4 STREAMING INSTABILITY 

Magnetic fields can be generated by streaming instab ility induced by cosmic ra ys, although 



alternative models have also been prop osed (see e .g. iGiacalone fc Jokipiil 120071) . Streaming 
instability can be induced in a resonant (lBelllll978l ) or non-resonant (lBellll2004l ) way, depend- 
ing on the type of interaction between particles and waves. In the former case, the unstable 
modes are Alfven waves, while in the latter case the modes are almost purely growing modes 
and do not correspond to Alfven waves. A satisfactory description of the interaction of par- 
ticles with these waves is missing, therefore it is very problematic at the present time to 
desc ribe cosmic ray diffusion in a background of wa ves which are excited non- resonantly 



see 



Zirakashvili fc Ptuskin 



2008a : 



Reville et al. 



2008 



for some recent attempts). Moreover, 



even the jump conditions at the subshock and in the precursor would be different from the 
ones typically adopt ed and it is not c lear a s yet which form the wave terms shuld have. In 



the context of SNRs, 



Amato fc Blasil (l2009l ) showed that the non-resonant modes are bound 



to be relevant only in the early stages of the evolution, while for most of the history of 
the SNR, streaming instability should be dominated by the excitation of resonant waves. 



For these reasons we chose to focus only on the self-generation of resonant A 



f ven w aves, 



leaving other cases, including the phenomenological approach of iBell Sz Lucekl (l200ll ). for 
future work. 

The stationary equa tion for the growth and transport of magnetic turbulence reads (e.g. 



McKenzie fc Volk 



19821 ): 
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uix] 



+ a{k, x)Vw{k, x) — T{k, x)Vw{k, x) 



(38) 



dx dx 

where Tw and Vw are, respectively, the energy flux and pressure per unit logarithmic band- 
width of waves with wavenumber /c, a{k,x) is the rate at which the energy in magnetic 
turbulence grows and T{k,x) is the rate at which it is damped. This equation is very gen- 
eral, but in the case of modified shocks one should keep in mind the fact that Fourier analysis 
in /c-space is only accurate for wavenumbers such that 1/k remains appreciably smaller than 
the typical length scale of the precursor. 

In the following we only consider resonant scattering between the accelerated particles 
and the magnetic turbulence, which gives for the growth-rate of energy in Alfven waves: 



a{k,x) 



47r va{x) 
3 Vw{k,x) 



p^v{p) 



dfix,p) 
dx 



(39) 



- p=p(fc) 



Amato &: Blasi 



2006, for a 



where p{k) = eB/knipC is the resonance condition (see e.g. 
derivation of this expression). 

Assuming no damping for the moment, and integrating Eq. [HE] in k-space we obtain: 



dF^{x) .dp^{x) dpc{x) 

= M(a;)— + va—: — . (40) 

dx dx dx 

We further assume that only waves of one sign of helicity are generated upstream and that 
Va ^ u, so that F^{x) ~ 3u{x)pw{x). With these assumptions Eq. HO] becomes: 



2U{x 



da{x) 



Va(x 



d^{x) 



3a(x 



dU{x) 



(41) 



dx dx dx 

Provided the cosmic ray generation is efficient and that both Mq and Ma are much larger 
than 1, we can neglect the plasma and the magnetic pressure with respect to the kinetic and 
cosmic ray terms in Eq. [351 hence ^(x) ~ 1 — U{x) and Eq. Hi] can be solved analytically, 
returning 

1 - U{xf' 



a{x) = f/(x)-='/' 



"0 + 



AMao 



(42) 



We could use Eq. UT] directly in our calculations, but it is easy to check that the assump- 
tions which lead to Eq. |42]are well satisfied in all cases of interest. The important physical 
information to retain at this level is the enhancement of the magnetic field due to adiabatic 
compression, which clearly shows in Eq. |42] through the factor U^^'"^. 

In the following we consider the case a^ = 0, assuming that all the turbulence is generated 
via streaming instability, thus 
1 - U{xf 



a{x) 



4Ma{x)U{x) 



(43) 
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Figure 1. Left panel: velocity (top) and cosmic ray pressure {bottom) profiles. Right panel: momentum distribution (top) and 
spectral slope (bottom) of the cosmic ray spectrum at the shock location. The thick lines represent the solution found when 
the magnetization is correctly treated in the jump conditions, while the thin lines show the solution found when it is ignored. 
Different line-types refer to different values of the Mach number Mq, as specified in the figure. The location upstream is in 
units of X* = —D{pmax)/uo and is expressed in base-10 logarithm, while the downstream is taken as homogeneous. 



where Ma{x) = Mao^/U{x) is the local Alfenic Mach number. 

For the sake of clarity, we stress that contributions of order 1/Ma are usually negligible, 
as in the calculation of the reflection and transmission coefficients or in the treatment of the 
magnetic energy flux. The only exception to this rule is in the conservation equations for 
momentum and energy, where the magnetic terms, of order 1/Ma, are comparable to the 
ones pertaining the gas, usually of order 1/Mq. Both these contributions are very relevant 
to the shock dynamics because they affect the compressibility of the system. However, they 
are extremely small compared to the kinetic and cosmic ray energy, so that Eq. |42] is a very 
good approximation. 



5 RESULTS 

In this section we show the results obtained through the algorithm described in Sec. El and 
including only resonant amplification of the magnetic field in the precursor, as described in 
SecH 

Here and in the following, unless specified otherwise, we assume a SNR age of lOOOyr, a 
circumstellar density of po = 0.5mp/cm^ and a background magnetic field of -Bq = 5|J.G. The 
injection parameter tp is kept fixed: ip = 3.7. We also consider two different circumstellar 
temperatures Tq = 10^ and 10^ K, which for mq = 5900km/s correspond to Mq = 500 
and Mq = 50 respectively. This is done in order to investigate different scenarios for shock 
propagation. 
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In Fig. [T]the velocity and cosmic ray pressure are shown {top and bottom panels on the 
left respectively) for Mq = 50 and 500 {dashed and solid lines respectively) with or without 
inclusion of the magnetic feedback {thick or thin lines). The values found for the relevant 
parameters in the same cases are reported in Tab. [2l 

The most striking effect that comes from the correct treatment of the jump conditions for 
the magnetic field is indeed the reduced shock modification, visible from the precursor profile: 
the total compression ratio is found to be Rfot ~ 9 for both values of the Mach number, 
while the prediction of standard non-linear theory would be Rtot ~ 112 for Mq = 500 and 
Rtot ~ 17 for Mo = 50. 

Clearly, the shallower precursor comes from the decrease of the fraction of bulk energy 
that is converted into cosmic rays, which however remains considerable: ^i is reduced from 
more than 90% to around 50-60%, for both values of Mq (Fig. [H left bottom panel). 

The remaining fraction of bulk energy ends up being converted into heating and energy 
of the turbulent magnetic field. We define the thermal emissivity as 

e{x) oc p{x)^T{x)^^^ oc U{xy^T{xf/^ , (44) 

and we report, in the last column of Tab. [21 the ratio between its downstream value com- 
puted within modified shock theory {52) and the prediction for the same shock in test particle 
theory {etp)- The temperature (next to last column in Tab. [2]) and thermal emissivity down- 
stream are always much smaller than in the test-particle approximation, since a considerable 
fraction of the shock ram pressure is now going into particle acceleration and magnetic field 
amplification. Nevertheless, it is clear that the effect of the magnetized jump conditions is 
to enhance both quantities (e.g. T2 is enhanced by a factor > 100 for Tq = 10^ K), as a net 
result of the increased temperature jump T2/T1 (Eqs. fTTfTOj) and the reduced compressibility 
of the upstream plasma (lower -Rtot)- 

We should recall that the compression ratios actually felt by cosmic rays depend on the 
relative velocity between them and the scattering centers, i.e. the Alfven waves, according to 
Eq. [22I Nevertheless, when the Alfven velocity is taken according to Eq. [211 the discrepancy 
between Sgub — Stot and Rgub — Rtot is usually negligible, the difference in the plasma's and 
cosmic rays' compression ratios being of order 1% at most for the cases in Tab. [2l 

As to the magnetic field, the values of B2 in Tab. [2] are in the correct range inferred from 
fits to X-ray observations of SNRs, which indicate B2 ~ 400 — 500|J.G for SNRs with uq as 
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Table 2. Comparison of the results obtained with and without inclusion of the magnetic feedback (As) for the 
shocks of Fig. [T] 



ro(K) 


As 


«i 


Pma:.(106GcV/c) 


Rsub 


Rtot 


^sub 


Stot 


Bii^G) 


TaClO^K) 


e2/etp 


10-* 
10* 


No 
Yes 


0.97 
0.58 


0.24 

1.17 


3.58 
3.84 


112.1 
9.22 


3.43 
3.79 


108.7 
9.12 


645.8 
463.9 


0.88 
126.5 


0.030 
0.346 


10*^ 
10"^ 


No 
Yes 


0.77 
0.54 


0.59 
1.14 


3.76 
3.84 


16.6 

8.44 


3.70 
3.79 


16.4 
8.36 


235.0 
425.1 


42.3 
154.8 


0.216 
0.391 



high as 5000 — 6000kni/s. This shows that amphfication due to streaming cosmic rays can 
actually account for such high magnetization levels. 

Let us now consider the spectrum of the accelerated particles for the same situations 
discussed above. These are plotted in the right top panel of Fig. [1] (distribution function in 
momentum space, multiplied by p^), while the right lower panel shows the spectral slope 
q{p) at the shock location. For diffusion coefficients that increase with particle momentum, 
higher energy particles stream further upstream of the shock than the less energetic ones. In 
the case of a modified shock, this causes the former to experience compression factors even 
larger than 4, the typical limit for strong shocks in the test particle regime. Since the spectral 
slope is basically determined by this effective compression ratio, the resulting spectra are 
concave, softer at the lowest energies and harder at the highest energies. It is intuitively 
clear that a smoothening of the precursor (i.e. a reduction of the ratio Rtot/ Rsub) results in a 
reduced concavity of the spectra as compared to the standard prediction of nonlinear models 
for a given Mach number of the shock. This effect is evident from the comparison of thick 
and thin curves in Fig. [H where the dynamical reaction of the amplified field is included: 
the particle spectrum is roughly close to p~^ up to lO'^GeV/c and tends to be fiatter above 
this energy, with a slowly changing slope. 

Apart from this global concavity of the spectra, three main differences between the case 
with and without the correct jump conditions are worth being noticed: 

1) Pinj increases due to the more efficient heating of the downstream plasma, and the 
fraction of particles injected into the acceleration mechanism (Eq. [25]) also slightly increases 
(from rj = 1.03 x 10""^ to r/ = 1.19 x 10""^ for Mq = 500). This is consistent with the increase 
in Rsub induced by the magnetic feedback. 

2) at the highest energies the spectra are somewhat steeper than usually predicted for 
strongly modified shocks (p^^-^ instead of p"^'"^), but the total cosmic ray energy is still 
dominated by the highest momenta; 

3) the maximum momentum achieved by non-thermal particles is increased by about a 
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Table 3. Solution of DSA for different SNR environmental parameters. 



To(K) 


no (cm 3) 


Boi^G) 


€i 


Pma^ilO^G 


=V/c) 


Rsub 


Rtot 


W 


Bii^G) 


T2ilO^K) 


10* 


1 


5 


0.70 


1.42 




3.71 


12.5 


141.9 


714.2 


61.3 


10* 


0.5 


10 


0.54 


1.60 




3.76 


8.18 


373.3 


579.0 


149.8 


10-* 


0.1 


5 


0.50 


0.74 




3.78 


7.68 


406.8 


255.1 


173.6 


10-* 


0.1 


1 


0.74 


0.36 




3.74 


14.3 


89.9 


201.9 


48.0 



10^ 


1 


5 


0.65 


1.43 


3.73 


10.8 


1.41 


632.2 


88.2 


10'^ 


0.5 


10 


0.51 


1.56 


3.77 


7.75 


3.64 


546.6 


171.5 


10^ 


0.1 


5 


0.48 


0.72 


3.78 


7.36 


4.00 


243.8 


192.9 


106 


0.1 


1 


0.67 


0.36 


3.78 


11.4 


0.89 


167.8 


82.5 


10^ 


0.01 


5 


0.12 


0.15 


3.94 


4.50 


4.41 


54.7 


586.0 


10^ 


0.01 


1 


0.37 


0.16 


3.92 


6.24 


2.17 


50.2 


302.6 



factor 2-5, and reaches a few times lO^GeV/c, as a result of the smoothening of the precursor. 
We point out that during the free expansion phase the shock velocity may be high enough 
to make the non resonant streaming instability substantially contribute to the magnetic 
field amplification. This would reflect in a boost of Pmax for two concurring reasons: first, 
as a direct consequence of the resulting decrease of the diffusion coefficient; second, as a 
consequence of an enhanced magnetic feedback that would cause a further reduction of Rtot 
(see Eq. [31] for the time needed to accelerate a particle to Pmax)- For a typical remnant age 
of r = lOOOyr, the maximum energies we derive appear to be in qualitative agreemen t with 



the knee in the proton spectrum as observed e.g. by KASCADE (jAntoni et al.ll2005l ). 

Finally, we notice that our results show a weak dependence on the background temper- 
ature To (in the range 10^ — 10^ K), and therefore on the sonic Mach number, as long as 
this is much larger than 1. This contrasts with the standard non linear prediction for the 
increase of the shock modification as Rtot oc Mq . The explanation lies in the fact that the 
magnetic backreaction is much more effective for high Mq, since it is driven by 

rs / D \ 7-3/2 



W - — --^^ / ^s«6 



Pi 4 Mao V R 



Hot 



p \ 2 



R, 



tot 






In Tab. [3] we show the results for several different choices of the environmental param- 
eters: it is clear that the magnetic feedback is never negligible (namely we always have 
ly > 1) and that the prediction for the compression ratios is consistent with the values 
inferred from observations. The case with Tq = 10^ K and low density no = O.Olcm"'^, which 
may be representative of the hot phase of the interstellar medium, does not lead to a shock 
as strongly modified as in the other cases because in such an environment an age older than 
~ lOOOyr is needed to achieve a strong modification. 
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5.1 The effect of the velocity of scattering centers on the spectrum 

As discussed in Sec. |3] , the propagation of the accelerated particles in the shock region is 
described by Eq. |23l The terms u{x) — va{x) describe the velocity of the scattering centers. 
Here va represents the wave velocity, which, as discussed above, we assume to be the Alfven 
speed calculated in the background field Bq. Provided the waves remain Alfven waves even 
in the regime 6B/B ^ 1 this result holds in an exact way. It is however clear that this 
can only be considered as an ideal case, in that turbulence may change such a picture in a 
considerable way. One way in which the change may appear is in changing the wave speed. 
In this section we investigate the effects on the spectrum of accelerated particles and on the 
shock precursor which derive from calculating the Alfven speed in the local amplified field, 
namely: 

..(.) = -^m.. (46) 

We stress that in our opinion this assumption is totally unjustified from the physical point 
of view, and we use this case only as a toy model to illustrate how sensitive the results can 
be to unknown non-linear effe cts. We notice however t hat sim ilar approaches have in fact 
been adopted , for i nstance by IZirakashvili fc PtuskinI (J2008bl ) and also, in some form, by 



BeU fc Lucekl fl200lh . 



The net effect of this apparently harmless assumption is that the velocity of the scattering 
centers is greatly enhanced and this affects in an substantial way the effective compression 



Belli fll978h 



factor at the subshock and therefore the spectrum, as was already pointed out by 
in the context of test particle theory. 

In order to illustrate this effect in a quantitative way we run our calculations for pi^j = 
3.7pi/i,2, Mq = 250, To = 10^ K, uq = O.Scm^'^, Bq = 5(J.G and an age of the system of 
~ lOOOyr. The results are illustrated in Fig. [2l The solid lines are obtained using the Alfven 
speed calculated in the background magnetic field Bq, as in the previous section. In this 
case the precursor is very evident (top left panel in Fig. [2]) and the spectrum of accelerated 
particles is concave (top right panel in Fig. [2]). 

The dash-dotted lines refer to the case in which the Alfven speed is calculated using the 
amplified magnetic field, but this definition is used only in the transport equation, while 
the rate of growth of unstable waves is left unchanged. As a consequence of the reduced 
effective compression ratio of the scattering centers' velocity felt by accelerated particles, 
the spectrum becomes softer, as illustrated in the right panels of Fig. [2l In fact one can 
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Figure 2. Left panel: velocity (top) and normalized magnetic energy density in the shock region (bottom) . Right panel: 
momentum distribution (top) and spectral slope (botttom) of the cosmic ray spectrum at the shock location. The solid lines 
refer to the solution found when the Alfven speed is calculated in the background field. The dash-dotted lines refer to the case 
in which the modified Alfven speed is used only in the transport equation, while the dashed lines refer to the case in which the 
modified Alfven speed is used both in the transport equation and in the growth rate of unstable modes. 



see that the spectrum becomes even steeper than p""^ at all momenta. When the modified 
Alfven speed is used also in the growth rate (dashed lines), the effect is even more dramatic 
and the spectrum has a slope ~ 4.3 at all momenta and the concavity is hardly visible. 
Moreover, the amplified field becomes smaller in the latter two cases, which also results in 
lower maximum energy of the accelerated particles, as visible in the right panels of Fig. [2j 
Despite the fact that the assumption of very large va could be unphysical, and certainly 
not well justified from the theoretical point of view, we cannot refrain from being very con- 
cerned by the dependence of the results on the value of the effective velocity of the scattering 
centers. On the other hand, having larger fields does not n ecessarily im ply f aster scattering 



Belli f l2004h and 



Amato fc Blasi 



center s. For instance, the non-resonant waves discussed by 

( I2OO9I ) are almost purely growing modes and one could expect that they may be almost 
stationary in the fiuid frame. In addition, one should keep in mind that if the velocity of 
the waves becomes too high, they may generate shocks in the background plasma and damp 
their energy on it, so that those waves do not contribute to the scattering of particles. 



6 THE (REDUCED) ROLE OF TURBULENT HEATING 



Herewe address another very important issue, already introduced in the work by iCaprioli et al 
( l2008al ). The correct inclusion of the turbulent magnetic field in the calculation of the jump 
conditions very naturally leads to values of Rtot that allow to fit the X-ray observations in a 
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totally new way in the perspective of cosmic ray modified shocks, i.e. without invoking the 
presence of additional gas heating mechanisms. 

It is well known that the shock dynamics, and in turn the particle acceleration process, is 
very sensitive to the presence of non-adiabatic gas heating, since a hotter upstream plasma 
is naturally less compressible. The low compression ratios inferred from observations have 
been often explain ed by invoking mechanisrns of non-adiabatic heati ng in the precursor, such 



as Alfven heating (jyolk fc McKenzie 



(JDrury fc Falle 



1986 



1981 



McKenzie fc Volk 



Wagner. Falle fc Hartquist 



2003). 



19821 ) and acoustic instability 



Acoustic instability develops when sound waves propagate in the pressure gradient in- 
duced by cosmic rays in the shock precursor. The instability results in the formation of weak 
shocks upstrear n, which cause heating in the prec ursor, thereby reducing the compressibility 
of the plasma. I Wagner. Falle fc HartquistI (120071 ) showed that there is however a range of 
steady state solutions characterized by moderate cosmic-ray acceleration and compression 
ratios significantly larger than 7. 

Alfven heating (also called turbulent heating) is a generic expression which is supposed, 
at least in principle, to apply to any damping mechanism for Alfven waves, and may result 
for instance from ion-neutral damping or non-linear Landau damping, depending on the 
ionization level of the ba ckground plasma. Although the initial mathematical approach of 



McKenzie fc Volkl (1l982l ) was based on the assumption of non-linear Landau damping, the 
formalism was generic enough that it could be adapted to any damping mechanism, and 
this is indeed what happened. The common formulation of the Alfven heating assumes that 
some fraction of the energy in the form of waves is damped into the thermal energy of 
the background plasma, independently of the details. Notice that this formalism does not 
distinguish among modes with different wavenumber k, therefore this type of calculations 
is intrinsically insensitive to the spectrum of turbulence and should not be used to infer 
information on the shape of the diffusion coefficient. Notice also that whenever applied 



to the case of non-linear Landau damping, t 
4000km/s(To/5 x lO^iTl^/^ Jvolk fc McKenzie 



le m echanism is effective only when Uq -C 



1981 



). It is not obvious that this condition 



holds in the young SNRs listed above, since typically Uq ~ 3000 — 6000km/s. 

A few other points are worth being mentioned: 1) Alfven heating was first introduced 
in order to avoid that the amplified magnetic field could reach the nonlinear regime. On 
the other hand the same formalism is now used even in situations in which 6B ^ B. Much 
care should be taken in using these expressions in the nonlinear regime. 2) The formalism 
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Table 4. Alfven heating effects 



c 


5i 


Pma AiO'^GeV/c) 


Rsub 


Rtot 


Bi/Bo 


W 


B2(nG) 


TiilO^K) 


£2/£tp 





0.60 


1.17 


3.76 


9.52 


25.3 


1.941 


475.6 


114.6 


0.317 


0.5 


0.66 


0.84 


3.65 


10.96 


20.8 


0.390 


379.6 


132.6 


0.523 


0.8 


0.65 


0.53 


3.68 


10.76 


12.8 


0.115 


232.5 


128.3 


0.480 


0.99 


0.55 


0.12 


3.85 


8.69 


2.26 


0.005 


43.5 


162.2 


0.553 



Solutions found for a shock with Tq = 10^ K, Mg = 50, po = Imp/cm^ and Bq = 5|a.G when the turbulent heating 
is taken into account according to Eq. 1501 for various ( < 1 (see Eq. I47I I . The last column shows the downstream 
thermal emissivity £2, normalized to the value it would have for the same shock in the test particle regime. 



introduced by lVolk fc McKenzid (Il98ll ) and iMcKenzie fc Volkl (Il982l ). and later adopted by 
virtually all authors willing to include Alfven heating, is based on the implicit assumption 
that there is a rapid damping of all energy of waves onto the background plasma (as stressed 
above, this was done in order to avoid excessive magnetic field amplication). This inhibits the 
generation of magnetic field, so that large magnetic fluctuations become unfeasible. In other 
words, the standard treatment of Alfven heating, which corresponds to having a{x) = T{x) 
in Eq. [3HI is incompatible with having large values of the turbulent field at the shock, as we 
show below in a quantitative way. 

In the following we assume that the damping rate is limited to a fraction of the growth 
rate, namely that 

r(a;) = Ca{x) . (47) 

An equation describing the Alfven heating of the precursor can be obtained by taking 



the derivative of the equation for the conservation of energy toget 
transport of waves and cosmic rays (see the derivation of Eq. 9 in 
Under the above assumptions we obtain: 

^ [PMxPixr] = c(7 - i)VAix)uixr-'^ , 



ler with the equations o f 



McKenzie fc Volk 



1982). 



(48) 



and Pth is now the pressure of the plasma including the effects of the turbulent heating. 
Clearly, in the limit of adiabatic evolution of the precursor, one has C = ^^^d Eq. [15] is 
recovered. 



It has been shown (JBerezhko fc Ellison 



1999 



Amato &: Blasi 



20061 ) that a good approx- 



imation to the solution of Eq. HH] is 

P^H(x)^P(x)|l + C(7-l)^[l-f/(x)l|. (49) 

where P{x) is the plasma pressure as calculated taking into account only adiabatic com- 
pression in the precursor, Eq. [151 This expression, which serves as an equation of state 
for the gas in the presence of effective Alfven heating, reduces to the standard Eq. 50 of 
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Berezhko fc EUisoru (119991 ) ioi C, = 1, while, for C < 1 the damping of the waves is mitigated 
and an effective amphfication of the magnetic field is allowed. 

The change in the equation of state of the gas also manifests itself in the shock dynamics. 



The new relation between the compression ratios Rsuh and Rtot reads 



R 



■7+1 
tot 



M'^R'^ 



7 + 1 - i?,„b(7 - 1) 
(1 + AB)il + Ath) 



where we have introduced 



A 



TH 



C(7-l) 



Ml 
Ma 



R 



sub 



R 



tot 



(50) 



(51) 



with a notation which allows a direct comparison between the effects of magnetic feedback 
and Alfven heating. 

It is widely known that the inclusion of Alfven heating has an important impact on 



the tota' 



[see e.g. 



3/4 



compression ratio , chang ing its scaling with the Mach number from Mq to M^ 



.3/8 



Berezhko fc Ellison 



19991 ). However the situation is very different when the correct 



magnetized jump conditions are taken into account. In Tab. HI we report, for different values 
of C = 0, 0.5, 0.8, 0.99, the solutions of the problem including the full treatment of growth 
and damping of Alfven waves according to the approximate analytical solution of Eq. [381 
along with the prescription of Eq. HH namely 

1-U{xf' 



axH 



(x) = (l-C) U{x)-^'^ 



"0 + 



AMaq 



(52) 



It is clear that the increasing relevance of Alfven heating {(^ approaching 1) does not lead, 
in this approach, to a smoother precursor (larger value of Rsub/Rtot)- In fact, the energy 
transfer from the waves to the plasma, while heating the plasma and reducing its compress- 
ibility, is also accompanied by a decrease of W, the ratio of magnetic/plasma pressure. The 
latter is the parameter which controls the magnetic feedback, which is then reduced. The 
net effect is a slight increase of Rtot for intermediate values of^ = 0.5 — 0.8. Only if ( is very 
close to 1 the shock is less modified than the adiabatic solution with magnetic backreaction, 
but even this effect is rather limited, since the decrease of Rtot is only about ~ 10% for the 
case C = 0.99. 

The main effect of the inclusion of the Alfven heating is instead a significant reduction 
of the magnetic field (more than a factor 10 between C = and ( = 0.99), which also leads 
to a correspondingly lower Pmax {Pmax — lO^GeV/c for ( = 0.99). 

We notice that the downstream temperature is affected in two different ways by the 
correct jump conditions and by the turbulent heating: the former provides an increase in 
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the jump T2/T1 proportional to W, while the latter results in a higher Ti. The interplay 
between the two effects produces the non-monotonic trend of T2 and of the thermal emissivity 
£2 (shown in the last columns of Tab. HI). The latter is slightly increased by the presence of 
non-adiabatic heating, with a value around 0.5etp, basically for all choices of C > 0. 



A recent investigation on the role o: 
modified shocks has been carried out by 



turbulent heating on the propertie s of cosmic ray 



Vladimirov. Ellison fc Bykovi (120081 ). These authors 



adopt a different recipe for particle injection into the acceleration process (which is a poorly 
understood issue in any case), that leads to a dependence of the fraction of injected particles 
on the temperature immediately upstream of the shock. As a consequence they find that 
dissipation of magnetic turbulence into heat upstream of the shock boosts the injection 
by a large factor and modifies the cosmic ray spectrum at low energies, although it does 
not significantly affect the overall acceleration efficiency and the high energy part of the 
spectrum unless the the fraction of turbulence energy that is transferred to heat gets close 
to 100%. If the latter is the case, on the contrary, also I Vladimirov. Ellison fc Bykovi (120081) 
observe a considerable decrease of both the magnetic field strength and Pmax, in agreement 
with our results. 

Before concluding this section, we summarize our main conclusions concerning the effects 
of turbulent heating: 



• if turbulent heating is in the form of nonlinear Land au damping, it is suppressed for 
uo > 4000km/s(To/5 x 10^ K)^/^ jvolk fc McKenzielll98l[ ). thus it could be important only 



if the circumstellar temperature is rather high (for instance because the shock is expanding 
in the pre-supernova stellar wind, or because of a dominance of the hot coronal gas phase) 
and the shock has already slowed down significantly; 

• if the rate of damping of the waves is too close to that of growth (C ~ 1), the magnetic 
field amplification is heavily suppressed and explaining the large levels of magnetization 
inferred from X-ray observations becomes very challenging; 

• if the damping is effective but still allowing sufficient magnetic field amplification, the 
smoothening of the precursor is roughly unchanged with respect to the results obtained in 
the case of magnetic feedback alone. 
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7 CONCLUSIONS 

The possibility that the narrow filaments detected in X-rays may result from severe syn- 
chrotron losses of high energy electrons has sparked much attention on the issue of magnetic 
field amplification at shock fronts, since the inferred fields are ~ 100 times larger than the 
typical interstellar ones. The importance of such magnetic fields for the origin of cosmic rays 
is immediate: if the turbulent field leads to Bohm diffusion, it is easy to show that particle 
acceleration at shock fronts may lead to the generation of protons with energy close to the 
knee. 



In a previous paper (ICaprioli et al 



2008al ) we investigated the hydrodynamical conse- 



quences of having a large magnetic field amplified by cosmic ray streaming upstream of 
the shock and we demonstrated that when the magnetic pressure equals or overcomes the 
pressure of the background plasma, a condition easy to realize, the compressibility of the 
plasma is reduced, so to decrease the compression factor of the cosmic ray modified shock 
and smoothen its precursor. 

In the present paper we went beyond the hydrodynamical picture and we solved the com- 
bined system of conservation equations, cosmic ray diffusion-convection equation, and equa- 
tion for the magnetic field amplification in order to investigate the effects of the precursor 
smoothening on the spectrum of accelerated particles. We find that resonant streaming in- 
stability is sufficient to amplify a parallel pre-existing magnetic field to levels which compare 
well with X-ray observations, if the latter are interpreted as a result of strong synchrotron 
losses. 

In these circumstances, we also confirm the crucial role of the dynamical feedback of 
the magnetic field, which leads to total compression factors around ~ 7 — 10 (to be com- 
pared with the typical predictions of standard NLDSA, which gives Rtot ~ 20 — 100), in 



good agreement with the va! 



ues suggested by observations and by fits to the mu 



spectrum of several SNRs (jVolk. Berezhko fc Ksenofontov 



2005 



Warren et a] 



tifreq uency 



20051). The 



reduced compression in the precursor does not inhibit particle acceleration, in fact a frac- 
tion 50-60% of the ram pressure at the shock is converted into accelerated particles, for 
parameters typical of a SNR in the Sedov phase. 

The kinetic calc ulation of partic l e accelerati on in the nonlinear regime was carried out us- 



ing the approach o: 
were treated as in 



Amato fc Blasil ( 



SkillineJ fll975bh 



2005, 



Bell 



2006). The effects of resonant streaming instability 



(119781 ) , but including the presence of the precursor 
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which iinphes magnetic field compression in addition__to amplification. We neglected here 



the possibility of having non- resonant amplification (JBelll l2004l ) . for several reasons: 1) the 



3), 



relation between pressure and energy fiux of these waves is not yet well defined; 2) the dif- 
fusion properties of charged particles in turbulence at wavelengths much shorter than the 
Larmor radius o f particles are no t know n as yet, and need a dedicated effort of investigation; 
3) as showed by lAmato fc Blasil (120091 ). the non resonant modes are likely to be especially 
important in the free expansion phase of a SNR and in the early stages of the Sedov phase. 

Perhaps not surprisingly, we find that the effect of magnetic feedback on the spectrum 
of accelerated particles is that of reducing the concavity. This result can be easily explained 
based on the reduced compressibility in the shock precursor, caused by the magnetic pressure. 
The slope of the spectrum typically remains close to 4 at energies smaller than ~ 1 TeV, 
while being flatter at higher energies (Fig. [1]). We also discuss the effect of the amplifled fleld 
on the velocity of the scattering centers, which causes a further steepening of the spectrum, 
leading it to get closer to a power law. We flnd that the spectrum of accelerated particles 
becomes even steeper than p~^ if the velocity of the scattering centers is assumed to be the 
Alfven speed in the amplifled fleld. Although this maybe an unphysical assumption, it serves 
the scope of showing that the results of calculations carried out in this strongly nonlinear 
regime may well be affected by details which at the present time we are unable to control. 

The smoother precursor induced by the magnetic feedback also implies a larger value of 
th e maximum achievable momentu m for the accelerated particles. We recall that, as showed 
by iBlasi. Amato &: Capriolil (120071 ) , a strong precursor leads to lower the maximum momen- 
tum since particles feel a smaller mean fluid speed and therefore a longer acceleration time. 



Using the acceleration time for modifled shocks as calculated by iBlasi. Amato &: Caprioli 
(120071 ). we flnd that for a lOOOyr old remnant the maximum energy can be as high as 
lO^GeV, a factor 5-10 larger than one could flnd without including the magnetic backreac- 
tion. In principle even larger maximum energies can be achieved by taking into account the 
non-resonant streaming instability in the early phases of the SNR evolution. 

Our entire analysis was carried out for the case of a shock propagating in the direc- 
tion parallel to the ambient magnetic fleld. The ampliflcation of the background fleld due 
to cosmic ray streaming soon changes this situation, in that the magnetic fleld upstream 
becomes highly turbulent. Changing the fleld obliquity should not change much the picture 
discussed in the present paper as long as the projection of the shock velocity along the 
fleld lines remains supersonic. If the supernova shock propagates in a medium in which the 
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background field is mainly oriented in a given direction, there will be locations at which the 
shock is quasi-parallel and others at which it is quasi-perpendicular. The implications for 
the morphology of the non-thermal emission require a dedicated investigation, in that the 
quasi-perpendicular regions could be expected to be bright because of the more efficient per- 
pendicular configuration, but quasi-parallel regions could be expected to be bright because 
of more efficient magnetic field amplification. Clearly the degree to which these expectations 
may be true depends on many details which at present are not under control, although 
observations of specific astrophysical objects (e.g. SN 1006) will provide us with important 
information on this issue. 

Finally, we investigated the effect of turbulent heating on the shock modification. We 



McKenzie fc Volk 



genera lized the formalism introduced by IVolk fc McKenzid (Il98ll ) and 
(1l982h in order to account for the possibility that a finite fraction of the wave energy is 
damped on the background plasma. We find that, although in general the heating of the 
upstream plasma results in a decrease of the compressibility, the effect is at most competitive, 
if not subdominant, compared with the dynamical feedback of the amplified magnetic field on 
the plasma. When turbulent heating is indeed dominant, a fraction very close to unity of the 
wave energy is transformed in thermal energy, and the main effect is that of suppressing the 
growth of the magnetic field, so that it becomes challenging to explain X-ray observations in 
terms of severe synchrotron losses. Moreover in this case the maximum energy of accelerated 
particles is drastically reduced and it becomes much smaller than the knee. 

It is worth stressing that while the magnetic feedback on the background plasma is based 
on well known physics, the turbulent heating is treated at best in a very phenomenological 
way, with little attention to the specific physical processes that may be responsible of the 
non-adiabatic heating and on how modes with different wavelengths are damped. In these 
circumstances we think that the effect of magnetic feedback is much more solidly assessed 
and should be considered as the chief mechanism for reducing the compression in cosmic ray 
modified shock waves with evidences for amplified magnetic fields. 

On the other hand, the role of turbulent heating may become much more important 
for older remnants, when the shock velocity drops below ~ 2000km/s. At these times, the 
amplification of the magnetic field may have become less important from the dynamical point 
of view, so that the magnetic feedback is also less important. These stages, as discussed by 



Ptuskin fc Zirakashvilil (120051 ) play an important role in determining the spectrum of cosmic 



rays observed at the Earth, at least at energies much lower than the knee energy. 
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